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P_l Abstract 

Two boundary value problems for an elliptic equation in divergence 

i-S^ form with bounded discontinuous coefficient are studied in a bidomain. 

On the interface, generalized dynamic boundary conditions such as of 

^ the Wentzell-type and Signorini-type transmission are considered in a 

subdifferential form. Several non-constant coefficients and nonlineari- 

7—i ties are the main objective of the present work. Generalized solutions 

^ are built via time discretization. 
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^ 1 Introduction 

• ^ 

^ In the description of real hfe phenomena, challenges in science and technology 

a such as diffusion problems with transmission conditions are being addressed 

(cf. for instance [7j and the references therein). We refer to [131 [E] a gen- 
eral framework which allows to prove, in a unified and systematic way, the 
analyticity of semigroups generated by operators with generalized Wentzell 
boundary conditions on function spaces with bounded trace operators. The 
thin obstacle problem (also called the Signorini problem) models thresh- 
old phenomena like contact problem, thermostatic device or semi-permeable 
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membranes [1]. In pQ the study relies on the presence of differential opera- 
tors. We point out that their method is based on a fixed point argument. 
Under continuous or even constant coefficients, the regularity was shown for 
the Laplace- Wentzell problem [12] or the thin obstacle problem 0. The 
question of dynamic boundary conditions can be found in frictional contact 
problems (see [20] and the references therein). Their theoretical and numer- 
ical achievements are based on the time discretization method being closely 
related to ours. 

With the aim of forcing to make realistic assumptions and then deal with 
the mathematical consequences, we prove the well-posedness of boundary- 
value problems subject to dynamic non-linear and friction-type boundary 
conditions. The present work extends the known results of Laplacian opera- 
tor to a general elliptic operator in divergence form with bounded measurable 
coefficient in the context of diffusion processes. The motivation comes es- 
sentially from the models for the electrical conduction in biological tissues 
[H El [TU] . The construction of generalized solutions is shown via time dis- 
cretization, following the Rothe method [lEl [TBI [19] . 

Let fli and be two disjoint bounded domains of ]R"(n > 2) such 
that = U O2 is connected with Lipschitz boundary. Let F = dQi fl 
Q C 80,2 denote a (n-l)-dimensional interface that can include the following 
descriptions. 

1. If dQi C Q then F is a closed curve (n=2) or surface {n > 3). Currently, 
Qi and ^2 are called the inner and the outer domains of Q, respectively. 

2. If Fl := Sr^i \ f = int(9fii f] dil) ^ ^ then 

• if n = 2, F is relatively open (see Fig. [l] (a)). 

• li n = 3, Qi stands for a cylindrical-type domain such that Fi 
represents its top and/or bottom (see Fig. [l] (b)). 

3. The case of dQi fl dQ ^ with meas(9i7i fl dQ) = can be clearly 
included whenever 9^2 is Lipschitz continuous (see Fig. [T](c)). 

In conclusion, we assume that dQi {i = 1,2) are Lipschitz continuous. The 
domains have neither cuts (cracks) nor cusps, and situations as in Fig. [l] 
(d) are excluded. Define a relatively open (n-l)-dimensional set F2 C 8^2 \ 
F, with meas(F2) > 0, and F/) = Fi U F2 where we will impose Dirichlet 
boundary conditions. 
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Figure 1: The geometry and interface conditions: 2D (a) and 3D (b) models 
when Fi ^ 0; (c) other possible situation; (d) 2D counterexample. 



Let us introduce the problems under study. For T > 0, find Ui : i7j x]0, T[— > 
M satisfying 

-V ■ (aiVui) = fi inn, {t = 1,2). (1) 

The first mathematical interest of this problem is due to the discontinuous 
coefficient which refiects the spatial dependence of the conductivity on the 
electrical conduction in different materials. 

On the exterior boundary dQ = {8^2 \ T) U Fi, we have homogeneous 
mixed boundary condition 



■ n = on dfl \ Yd and = on Yd. 



(2) 



On the interface F, we study two different types of dynamic bilateral condi- 
tions. 
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Wentzell-type transmission The generalized Wentzell transmission bound- 
ary condition is given by 

ui = U2 and (3) 
[aVu ■ n] + (3Aui - adtUi G dj{ui) on S := rx]0, T[, (4) 

under the initial condition 

Mi(-,0) = 5onr (5) 

where a and S are known functions and /3 is a non- negative constant. 
If /3 = 0, the transmission boundary condition (|3|-(|5]) looks for the 
transmission in a thin (or lower dimensional) porous layer. Here n is 
the normal unit vector to F pointing into Q2, d is the sub differential 
with respect to the argument of the function j, and [■] denotes the jump 
of a quantity across the interface in direction of n, e.g. [ctVm ■ n] := 
(T2V'U2 ■ n — ctiVmi ■ n. 

Signorini-type transmission The transmission that characterizes the bound- 
ary thin obstacle problems such as the semi-permeable membrane is 
constituted by the jump condition 

[aVu ■ n] = g on T, (6) 

and the Signorini-type boundary condition 

asVna ■ n - adt[u] E dj{[u]) on S := rx]0,T[, (7) 

accomplished with the initial condition 

[u](-,0) = 5 on r (8) 

where g, a, j and 5* are known functions [1]. 

The most common application is when dj represents the indicatrice Heav- 
iside. These boundary-value problems also model some of the slip phenomena 
observed in contact problems [HI |20] . Other related problems are the uni- 
lateral problems [3]. 

The paper is organized as follows. Next Section we set the functional 
space framework, the assumptions on the data and main results. Sections [3] 
and |6] are devoted to the proofs of existence and uniqueness of weak solutions 
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of each problem, namely provided by the Wentzell-type and Signorini-type 
transmission, respectively. These two Sections have similar structures based 
on the time-discretization technique and are split into several subsections 
in order to clarify the exposition. In Section [5} we show how the unique 
solution to the boundary value problem provided by a thin porous layer can 
be obtained as the limit of perturbed problems. Finally, some additional 
regularity is shown in corresponding Sections |4] and [7} 

2 Functional space framework and main re- 
sults 

The data are given under the following regularity assumptions. Here we 
assume that 

(Tj e : 3o-#,cr*>0, o-# < o-j(x) < o"*, for a.a.a; G fi^; (9) 

for z = 1,2, 

aeL°°{T): 3a#,a*>0, a#<a{x)<a*, for a.a.a; G T; (10) 
and j : M — 7- M is a convex and lower semicontinuous function such that 

j > and j(0) = 0. (11) 

Let us define 

HhM = {veH\n): v\r^ = 0}; 

Hl^{n,) = {veH\n,): v\r,=0}, (2 = 1,2). 

For a Lipschitz domain fii, the trace operator H^^{^li) — )■ Hqq'^{T) has 
bounded linear right inverse, that is, for every element S of the trace space 

iJoo^(r) = {v & -^^(r) : its zero extension belongs to H^^'^{dQi)} 

there exists G H^-^^Qi) such that = S on T [15]. However, the trace 
mapping considered as a mapping from H^^{Q2) in L'^{dVt2) is surjective on 

Hl!,\dn,\f,). 
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Considering that the Poincare inequahty occurs when (1 dQi ^ 0, for 
i = 1, 2, then the above Hilbert spaces are endowed with the norms 

Mm (Q,) = l|Vt;||2A- 

When Fi = and then we endow ifp^(fii) with any of the equivalent 
norms 

\\v\\2,ni + ||Vt;||2,f7i ~ ll^l|2,r + ||Vt;||2,ni- 
Indeed, we recognize that H^^i^i) = H\ni) and HQ^'^iT) = H^I'^{dVli). 

2.1 Wentzell-type transmission 

We can interpret the solutions Ui : i7jX]0,T[— )■ M (i = 1,2) as the uniquely 
(almost everywhere) determined function u : r2x]0,T[— )■ M such that u\^^ = 
Ml, u\n2 = U2 and ui = U2 on F. 

Let us define Hp as the Hilbert space 

{v G Hl^{n) : vi = v\n^; V2 = v\n^; Vi = V2 on F} if (3 = 0; 
{v e Hl^{n) : vi=v\n^; V2=v\n^] Vi=V2 on T; Vv e L\T)} if /3 > 0, 

endowed with the inner product 



{u, v)i3 = / Vn ■ Vvdx + l3 Vn ■ Vvds. 
Jq Jr 

Definition 2.1. We say that a function u G L'^{0,T; Hp) is a weak solution 
to the problem if dtU G i^^(S) and it satisfies ^ and the variational 

formulation 

(jVu ■ V{v - u)dxdt + (3 [ f Vu- V{v - u)dsdt + 



Jo. Jo Jr 



I adtu{v — u)dsdt + / / {j{v) — j{u)}dsdt > 
Jr Jo Jr 

> [ {f,v-u)ndt, yveL\^,T-Hp), (12) 
Jo 



with 
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The symbol (■, ■)q denotes the duahty pairing (-, ■)[Hfi)'xHfi- 
For u : fix]0,T[— > M such that the homogeneous Neumann boundary 
condition in ^ is satisfied, the Green formula yields 

-{W ■ {aWu),v)n= I ctVm ■ Vwdx + ([ctVm ■ n],t;)r, \/v e Hp. 
Jn 

Thus, using Q and Q it follows Q. 



Theorem 2.1. Under the assumptions 



eH(s: u° = S onT; (13) 

^j('5M^<C^(l|5||lr + l), (14) 

where C stands for a positive constant, and f G C^'^{0,T; {Hp)') with the 
Lipschitz constant d, that is, 

\\f{r) - /(t)ll(H,)' < d\T - tl Vr, t G]0, r[, (15) 

there exists u G L^{0,T; Hp) a unique weak solution in accordance to Defi- 
nition \2.1\ 



Remark 2.1. The assumption (I4) yields if for instance j verifies j{d) < 



C{d'^ + 1) for all d G M. Notice that (13) guarantees that S G -^^^(r) is such 
that I3WS G L2(r). 



Theorem 2.2. Let the assumptions of Theorem 2J_ be fulfilled. Moreover, 
if the compatibility condition 

[ aVu^ ■ V{v - u^)dx + (3 f Vu^ ■ V{v - u^)ds + I {]{v) - ]{S)}ds > 
Jn Jt Jr 

> {f(0),v~u')n (16) 

holds for ally G Hp, then dtu G L'^{0,T; Hp)nL^{0,T; L'^{T)). In particular, 
ueC{[0,T];Hp). 

The transmission problem in a thin porous layer, Q-Q with (3 = 0, can 
be obtained as the asymptotic limit, when a small parameter e goes to zero, 
of the following perturbed problem, whenever the interface F = dQi C Q, 
Fi = and F^ = 
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(Pg) Find Me : = f2i U S'e U Q2,e 1^ satisfying 

- (TiAue = /i in Qi] 
-cTsAue = /a in ^s.e; 
E'jAus — adtUe G dj{ue) in S'eX]0,r[; 

M,(-,0) = M° in S",; (17) 
[Ue] = [cVwe ■ n] = on F; 
K] = [aVw, • n] = on F, := \ F; 
Vm2 ■ n = on 9fi \ Fg; 
U2 = on F2, 

with Se = {^ + rn(0 : ^ G F, < r < £7(0} where 7_G C°'i(F) such that 
< 7# < 7(0 < 7* for all ^ G F, and e > such that S~, C n. 
Let us define the Hilbert space 

Xe = {v e H^^{ne) : vi = v\n^, vs, = v\s,, V2,e = v\n2.,; 

vi = vs, on F, vs, = V2,e on FJ, 

where f^^ = f2i U S'^ U ^2,e- 

Proposition 2.1. Let the assumptions ^-(11), |7g| ), (15) and (3 = be 
fulfilled, and {I4) be replaced by j{d) < C((i^ + 1) for all d G M. Then 



the unique solution u of the problem ([ijj-([5p in accordance to Theorem 2.1 
under the admissible test function space X := L'^{0,T; H) H H^{0,T; H^{Q\ 
is the limit of the sequence of the unique solutions to the variational 
formulation of the perturbed problem (P^) 

/ / aeVue -Viv — Ue)dxdt + i i — dtUs{v — Us)dxdt + 
Jo Jn^ Jo Js^ ^7 

+ / / —{3{v)-3{ue)}dxdt>[ {fe,v-u,)nJt,yveL\0,T;X,),{18) 
Jo Js, n Jo 

with {7^, = aixn, + Xs^ + (^2X^2,. o-'^d f^ = fiXn^ + f2Xn2,e- 
2.2 Signorini-type transmission 

Here, we keep the notation of jump [v] = f 2 — f 1 for any vector v = (f 1, f2). 
However, in order to differentiate this case from the above, let us set every 
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vector by boldface. In general if Vi 7^ V2 on F, their weak derivatives do not 
exist. Let us define the Hilbert space 



endowed with the norm (cf. Lemma 6.1 ) 

||v||v = ||Vi;i||2,Qi + ||Vt;2||2,Q2 + IIHIh.r- 
For V G V, v|r G i/o'o'(F) x hH^OQ^ \ Ts). 

Definition 2.2. We say that a function u = {ui,U2) G L^(0,T; V) is a weak 
solution to the problem (^-(^ with @-(^ if dt[u] G i^^(S) and it satisfies 
^ and the variational formulation 

r 

aVu • V(v — u)dxdt + / {g, Vi — Ui)rdt + 
'0 Jn Jo 

+ 1 I adt[u]i[v] - [u])dsdt + I l{j{[v])-j{[u])}dsdt> 
Jo Jr Jo Jr 

> [ (f,v-u)nrft, VvGL2(0,T;V), (19) 

with 

o- = (TiXni + (^2Xn2 and f = (/i, /2)- 

Here, we use the same notation {■,-)n to denote the duality pairing 
(■,-)v'xv- The symbol (■,-)r stands for the duality pairing {■,-)y'xy, using 
the notation Y = Hqq^IT). 

For u = {ui,U2) such that the homogeneous Neumann boundary condi- 
tion in ([2]) is satisfied, the Green formula yields 

-(V ■ (aVu),v)n = / aVu-Vvdx + {[aVu-n],Vi)r + {a2Vu2-n,[v])r, 
Jn 

for all V G V. Thus, using Q and Q-Q it follows (|T9[. 

Theorem 2.3. Assuming (^-([77p, p4^ , f and g are Lipschitz functions in 
the following sense: there exist two positive constants di and d2 such that 

||f(T) -f(t)||v' < (20) 
\\g{r)-g{t)\\Y' <d2\T-tl Vr,tG]0,r[. (21) 
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and 

3u° G V : = S onT, (22) 

there exists u G L°°(0,T; V) a unique weak solution in accordance to Defini- 
tionlKM 



Remark 2.2. The assumption (22) implies that 



||^||2,r< ||K]||2,r< ||u°||v. 



Theorem 2.4. Let the assumptions of Theorem 2^ be fulfilled. Moreover, 
if the compatibility condition 

aVu°-V(v-u°)cix+(^7(0),t;i-<)r+ J {ji[v])-jiS)}ds > (f (0), v-uO)^,, 

(23) 

holds for all \ e Y , then dtu e L^{0,T;Y) n L°°{0,T; L'^(T)) . In particular, 
ugC([0,T];V). 



3 Proof of Theorem 2.1 



3.1 Discretization in time 

In the foUowing we use similar arguments from the methods described in 
|18j . We decompose the time interval / = [0, T] into m subintervals li^m = 
[ti^m,ti+i,m] of size h = T/m, i G {0, 1, ■ ■ •, m — 1}, m G N. We define, for 
all i e {0, 1, ■ ■ m - 1}, = u{ti+ i^rn) as solutions given at the following 
Proposition. 

Proposition 3.1. Let z G {0, 1, ■ ■ -, m — 1} be fixed, G L'^{T), and 

= /(t.+i,™) G (Hp)'. 
Then there exists m'"*"^ G Hjj a solution to the problem 

[ aVu'^^ ■ Viv - u'+^)dx + (3 [ Vm*+^ ■ V(t; - u'+^)ds + 
Jn Jr 

+ ^u'+^ {v - u'+^) ds + j{j{v) - j{u'+^)]ds > 
> {f'+\ V - u'+')n + ^ ^u' {v - u'+') ds, Wv G H^. (24) 
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Proof. The existence of a solution to ( 24 ) is deduced from the general theory 



on maximal monotone mappings applied to elliptic variational inequalities 
[211 pp. 874-875, 892-893]. Indeed, the mapping A:Hp^ {Hp)' defined by 



is single-valued, linear and hemicontinuous; the mapping ip : Hp — )■ [0, +oo] 
defined by 

j^j{v)ds, iij{v)eL\T) 
+00, otherwise 



is convex, lower semicontinuous and ip ^ +oo; and the coercivity condition 
is valid 

{Au,u) + (p{u) = I a\Vu'^dx + (5 I |Vup(is > min{cr#, , 
under the assumptions ([9])-(ll). Then, for b E (Hp)' such that 



the variational inequality (24) has a unique weak solution u = G Ha. □ 



Remark 3.1. Since u = S onT means thatu G L iV), then Proposition 3.1 
guarantees the existence of E V and consequently G -Z^^(r). Therefore, 
Proposition 3. 1 successively guarantees the existence of u^^^ G V for every 



i = 1, • • • , m — 1. 

3.2 Existence of a limit u 

Proposition 3.2. For all i G {0, 1, ■ ■ ■ ,m — 1}, the estimate holds 

a#\\u'^%r < max{-^, l}||/||i.(o,T;(/.,)') + «*ll^ll2,r- (25) 

Moreover, if {um}meN is the sequence defined by the step functions 
Hp 

~ _ j ""^ /^'^ t = 

Um[t) — S i+i . -i, , 1 

then there exists u such that 

Um ^ u in L^{0, T; Hp). 
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Proof. Choosing f = as a test function in ( 24 ) , we get 
for alH G {0, 1, ■ ■ ■, m — 1}. Then it follows 

min{a#, 1}\\u^^Th, + l{u^^'fds < max{^, l}||f +1?^,), + ^ ^(n^^d.. 
Summing on k = 0, it follows 



k=0 k=l 



Consequently, we get (25) and, for i = m — 1, 

min{cr#, l}||Mm||i2(o,T;/f^) < max{-^, l}||/l|i2(o,T;(H^)') + «*il'^ll2,r- (26) 

Thus we can extract a subsequence, still denoted by Um, weakly convergent 
to u e L'^{0,T; His). □ 

Next, let us study the discrete derivative with respect to t at the time 

^ -~h — 
Proposition 3.3. Let Zm '■ [0,T[— > L^(fi) he defined by 

( fort = 

Zm{t) = \ in n. 

^ Z in ]ti,m; '^j+l,m] 

If the assumptions ^-(11) and (13)- (15) are fulfilled, then the estimate holds 



\'^m\\L°°{0,T;Hfi) + ll^mL.S ^ ^ ill f W L'2{0,T;{Hpy) + 11^ lli?^)" (^7) 



Hence, we can extract a subsequence, still denoted by Z^, weakly convergent 
to Z e 
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Proof. For a fixed t, there exists i G {0, ■ ■ ■, m — 1} such that t E.]ti^rn] ^i+i,m] 
Choosing f = as a test function in (24), we have 



a 



n Jt 



In order to sum the above expression on k = 0,...,i, consider the relation 
2{a — b)a = + (a — 6)^ — 6^ to obtain 



J2 f (^Vm*^+^ ■ V(7 



V / Vn^+i-V(n'=+^-M'=)cis = ^ / |VM*+iprfs- J / \Vu'\^ds + 
k=o -^r 2 Jr 2 

Now, using the assumptions ([9])-(ll) we find 
min{a#, 1} 



k=0 



•^^ k=l 



«'+0n.(28) 



By (15) it follows 

i i 



k=l 



k=l 



Therefore, inserting the above inequality in (28) and applying (26), it results 
(|27l). □ 
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From the Rothe function defined by 

uiix,t) = u\x)+t ^ ^ , — m Joi = /, 

h 

consider the following definition. 

Definition 3.1. We say that {um}mm is the Rothe sequence if 
for a// i G {0, 1, ■ ■ ■ , m — 1}. 



Proposition 3.4. If Z satisfies Proposition^^ then 

dfU = Z in L'^(T), for almost all t E I. 

Proof. For a fixed t, there exists i G {0, ■ ■ ■, m — 1} such that t G]tj,m; ii+i,m\- 
Thus we obtain 

I Z„i{T)dT = 2^ / r + / r '^'T ^• 

Jo Jkh h Jih h 

Because there exists w G C([0, T]; L'^{T)) such that 

{w{t),v)= [\ziT),v)dT, 

Jo 



3.1 



on r. Thus we have j^Zm{T)dT = Um{t) — S 



let us consider Definition 
and from the Riesz theorem we get 

{um{t)-S,v)= I {Zm{r),v)dT, \/veL\T). 



Indeed, the right hand side of the above equation is a bounded linear func- 
tional in i^^(r), representable thus (uniquely) by the element Umit) — S from 
L\T). 

Then it follows 

lim {u^it) -S - w{t),v) = lim / (Z„(r) - Z(r), v)dT = 0. (29) 

14 



Let us prove that the norms of the functions uniformly bounded 



with respect to t G / and m. From the estimates (25) independent on i and 
m, and considering 

II /i\ii II? 11,^ ^i.m \ , j' j_i ^ ti,m n 

\\uUt)hr = \\u' (1 + — j^j — ^"^'^ 

then, we get 

<C(II/II + ll^lllr)- 

Hence, the Lebesgue Dominated Convergence Theorem can be apphed in 



(29) giving 

T 



hm / (umit)- S -w{t),v)dt = 0, \/v e L\T). 

m— >+oo 







In the same manner this result can be derived for the case when v{t) is 
a piecewise constant function of t G /. Since these functions are dense in 
it remains valid for every function v G From the uniqueness of 

the weak limit, we conclude 

U{t)-S= [ Z{T)dT, 

Jo 

which corresponds to the claim. □ 
3.3 Passage to the limit on m -\-oo 

Denoting /^(t) = /*+^ for t e]ti^rn,ti+i,m] and i G {0, ■ ■ - 1}, we have 



crVum ■ Vvdxdt + (3 Vwrn ■ Vvdsdt + / aZmvdsdt + 
Jt, Jt. 

+ / j{v)dsdt> / a\Vum\'^dxdt + j3 / \Vum\'^dsdt + 
Jt. Jq Js 

+ / aZmUmdsdt+ / j{um)dsdt+ / {fm,v - Um)ndt. 



From Propositions |3.2| and |3.3| to pass to the limit the above inequality 
and recalling the weak lower semicontinuity property for the first and second 
terms on the right hand side of the above inequality, it remains to prove that 



Mm — M in L^(E). 
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Taking 



u = u 



m ""m 



u first let us prove that 



Mm — — ^ in L (S). 
Since we have < t — tj^m < h m. ti+i,m\ "we obtain 

\\Um{t) - Um{t)\\2,r = \\Zm\\2,r{,h - {t - U^ra)) < h\\Zm\\2,T 



and from (27) then it follows 

CT 



\Ur. 



Mm 2,E < 



m 



i^(o,w) + ll«lH,)'/'->o. 



Secondly the Rothe sequence {um\ is bounded in L^(0, T; if^), and, from 
Proposition 3.4, the functions dtUm are bounded in -Z^^(S) then, for a subse- 
quence still denoted by Um, the strong convergence holds 

Um — > u in L^(S). 

Then it results 



ZmUmdsdt — 7- 



Uudsdt 



dtuudsdt. 



Therefore we are in the conditions to pass to the limit concluding the 
weak formulation ( 12 ) . 

From the standard technique to prove uniqueness of solution, the solution 
u to (12) with (Is]) is unique. Then the whole sequence {um} converges *- 



weakly to u G L'^{0,T;H,3 



4 Regularity in time 



Proof of Theorem |2.2| . The proof follows the time discretization argu- 
ment as in Theorem |2.1 , considering the existence of the integral inequality 



(24). Choosing v 



test function in rt24| for the solutions 



u 



and M*, summing the consecutive integral inequalities, and dividing by 



h, we deduce 

ha\VZ'+^\'^dx + hf3 [ |VZ^+i|2 



ds 



a{Z'+^ - Z')Z'+^ds < 
<{r+'-f\Z'+')n 
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taking the convexity of j into account. Applying the assumptions ^ and 
(15), it results 



min{a#, ||L + / - Z')Z'+'ds < dh\\Z 



i+l 'yi\ ryi^X 



71+1 1 



Considering the relation 2 (a — 6) a = a"^ + {a — b)"^ — b"^ , to a = and 
b = Z\ and summing on /c = 1, ■ ■ ■ , z (i G {1, • ■ • , m — 1}) we obtain 



mm{a#, 1} h\\Z'^'\\H, + c^#\\Z'^%r < 2 / 
k=i -^r 



a 



h 



ds + 



1 * 
max{ — , 1} /i. 

fc=0 



Notice that mh = T. 

Let us determine the estimate for the first term on the right hand side 
of the above inequality. Rewrite the integral inequality ( 24 ) for i = in the 
form 

/ (tV(m^ - M°) ■ V{v - u^)dx + / o-Vm° ■ V{v - u^)dx + 
Jn Jq 

+(3 j V(m^ - M°) ■ V(f - u^)ds + /3 j Vu"- V{v - u^)ds + 
+ I a'^^{v-u')ds+ [{j(v)-j{u')}ds>{f'-f{0),v-u% + 



h 



for all V and in particular v = vP . Thus, we apply the assumption (16) 
with V = and divide by h we deduce 



2h 



Viv} -v^)\^dx + 



2h 



Viu' -u'')\'ds+ I a 
C 



h 



ds < 



<^\\f-fmUy 



Then, using (15), we have 



a 



h 



ds < Chd^ < C. 
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Since the above regularity estimates are independent on m the proof of 
the passage to the hmit is similar to the one of Section [3j Moreover, the 
uniqueness of the weak solution implies that the weak solution is the strong 
solution in the sense u G C{[0,T]; Hjs) by appealing to the Aubin-Lions 
Theorem. 



5 Proof of Proposition |2.1 

5.1 Existence of Up 



The time discretization described in Section |3.1| reads, for the perturbed 
problem, as 

[ aeVu'^' ■V(v-u'^')dx+ I -^{u'+^ - u') iv - u'+^) dx + 

+ [ -{3iv)-jiu^+')}dx>{r+\v-u^+%^, yveX,. (30) 

The existence and uniqueness of a solution ■u*^^ = G is due to 
standard results for elliptic variational inequalities as in the proof of Propo- 
sition 3.1| (cf. |17]). Indeed, the bilinear symmetric form 

a(u, v) = a^Vu ■ Vvdx + / — — uvdx 

is coercive in the following sense 

a{u,u) > min{l,(T#}||Vn||^^f^^ + -^^Mls,- 



Now taking first v = in (30), analogously to the proof of Proposition 



3.2, we get the estimates 



r ^ a* 1 

min{a#,l} / \\um\\\dt < - — +max{ — A}\\fe\\l2(^a,Tix,)'y (31) 
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Next taking f = wMn (30) and arguing as the proof of Proposition 
we obtain 



3.3 



min{l,a#}||Vn^^^||^,,^ + ^^||Z'=+^||^,,^< 



fc=0 



< 



-^jiu')dx + C{\\Vum^^ + ||/.||i.(o,T:(x.y) + -\\A\lsJ- 



Thus applying (14) it resuhs that and are uniformly bounded in 
L°°(0,T;Xe) and L^(5'£x]0, T[), respectively. Therefore the existence of a 
solution u G lv^(0, T; X^) to (18) can be done by similar arguments of passage 
to the limit as in the proof of Theorem 2.1 (cf. Section 3.3). 



5.2 Passage to the limit on e 

In order to let e — )■ 0, we utilize the following equivalent variational inequal- 
ities to (18) and (12) with /3 = 0, respectively, 

/ / aeVus -V {v — Ue)dxdt + i i — dtv{v — Ue)dxdt + 
Jo Jn^ Jo J ^1 



+ 



a 

2e^ 



\v{fi) -uydx + 



T 



Js, ^7 



-j{ue)}dxdt > 



> I {fe,v-u,)njt, \/veXe:=L^iO,T;Xe)nH\0,T;H\Se)); (32) 
Jo 



and 



Jn 



aVu ■ V{v — u)dxdt + 







a , 



adtv{v — u)dsdt + j (0) — + 



Jr 



+ / / {j(^) - j{u)}dsdt > {f,v- u)ndt, \/v e X. 



Let Ue be the solution of (18), or equivalently (32), satisfying (17). By 
appealing to Section we have 

ll^e||L°°{0,T;L2(S£)) < C*! || ""^ || 2,f7 + || / 1| L2(0,T;H') ) ■ 

Using the result (cf. [8J) 

hur.,s.<c{\\ur.,v+e\\yur.,s.) 
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in the estimate (31) it follows 

\\u£\\L^O,T;H^ne)) < ||/||l2(o,T;H'))- 

Thus there exists a subsequence £ — )■ and a function u E L°°{0, T; L^(5'e))n 
L^{0,T;H^{n^)) such that 

Ue^u *-weakly in L°°(0, T; L^Se)); (33) 
^ M weakly in L^{0,T; H^Qs))- (34) 

Next we recall the following lemma which is an extension the one proved 
in [a [9]. 

Lemma 5.1. a) For any function w G iy"'^'^(r2 \ Qi) we have 

w 

— dx — )■ / wds as £ 0. 



b) For any sequence of functions G L^{(yL \ r2i)x]0,T[) and any w G 
L^iTx]0,T[) such that 

rT 



|Vwe||q,s. < C and 



[Wi; — w)dsdt — )■ 0, 



for some constant C > and some exponent g > 1, we have 



dxdt 



Js, ^7 



wdsdt 



as e 



5.1 



For an arbitrary v e Xr ^ X^n C([0, T]]H^{Vt \ Vti)), by Lemma 
we have 

/ —\v{Q)-uydx^ I -\v{0) - u^l^ds. 
Js, 2^7 Jr 2 

In order to apply Lemma 5.1 b), we define Wg = {v — Ue)dtv and w = 

(f — u)dtv. By (33) we obtain 

/ / — w)dsdt — )• 0. 
Jo Jr 

Since dtVv G L^{Qx]0,T[) we have 

||Viye||g,s, < W'^iv - Ue)\\2,sJ\dtV\\^ s+ \\v - UsW^ sM^t'^^hSe 

2— q'^ 2— g'^ 

for g > 1 satisfying 2g/(2 — g) <2n/{n ~ 2) that means g < n/(n — 1). 

Thus we can pass to the limit on £ —!■ in (32) to obtain the desired 
solution. 
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6 Proof of Theorem 12.3 



The generalized version of the Poincare inequahty apphed to functions ad- 
mitting jumps [2] can once more extended to the following version. 

Lemma 6.1. Let v G V. Then 

[ vldx <C\ f |Vv|2rfx + f [vfdsX . (35) 



Proof. If Fi 7^ 0, the classical Poincare inequality is valid and then (35) 
clearly holds. If Fi = we will prove by contradiction. Assuming that (35) 
is not true, there exists a sequence {vm} C V such that for all m G N 

IK'imlb.Hi = 1 and ||Vv„|||f^ + IKt'JUir < ^M- 

Hence Vv^, — )■ in L^(fi) and [fm] — )■ in L'^iT)- Since V is a reflexive 
Banach space, we can extract a subsequence of v^, still denoted by v^, such 
that Vm ^ V in V. Thus Vv = in f2 and vi = V2 on F. Consequently 
Vi G ifp^(i7i) and V2 G //^^(^s) satisfy t>i = t>2 = 0. From the compact 
embedding V L'^iyii) x L'^{yi2) it follows that 

^ in L'^iVLi) x ^^(fia). 

Then we conclude that 

il^^lmlb,^! = 1 ||0||2,f7i = 1 

which is a contradiction. □ 



6.1 Discretization in time 

we will construct weak solutions u*'^^ 



3.1 



As in Section 

{0, 1, ■ ■ ■, m — 1}, of an approximate time-discrete problem. 



I G 



Proposition 6.1. Let the assumptions ^-{11) he valid, m > a^^T/a^^ and 
i G {0, 1, ■ ■ -, m - 1} he fixed, [u'] G ^^(F), 



pi+l 



f(ti+i,m) e V and g 



i+l 



9{t 



eY'. 
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Then there exists a time-discrete solution u*+^ G V to the problem 

[ aVu'+^ • V(v - u'+^)dx + / (H - + 

Jn Jr 



+ ll[^^']{[v]-[u'-''])ds, VvGV, (36) 



with [u^] = S on r. 



Proof. We show the existence of a solution to (36) with the aid of the gen- 
eral theory on maximal monotone mappings applied to elliptic variational 
inequalities |[21:, pp. 874-875, 892-893]. To this end, we define the mapping 
v4 : V ^ V by 

(y4u,v)= / aVu-Vvdx + / -MMrfs 
Jn Jr h 



which is single- valued, linear and hemicontinuous; and the mapping 99 : V ^ 

[0, +00] by 

+00, otherwise 



which is convex, lower semicontinuous and ^ +00. Because of ([9])-(ll) the 
coercivity condition 

(Au, u) + y9(u) = / a"|Vu|^(ix+ / —[ufds+ i j{[u])ds > a^\\\i\\'^, 
Jn Jr Jr 

is valid for any h < a^/a^. Then, for b G V such that 



(b, v) = -(r+\ v)f, + {g'+\v,)r - I -[u'][v]ds 



a , 



the variational inequality (36) has a unique weak solution u = u'"*"^ G V. □ 

Remark 6.1. Since [u^] = S on T means that [u^] G -^^^(r), then Proposi- 
tion 6.1 guarantees the existence of G V and consequently [u^] G -Z^^(r). 
Therefore, Proposition 6A_ successively guarantees the existence of vl'^^ G V 
for every 2 = 1, ■■■ ,m — 1. 
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6.2 Existence of a limit u 



Proposition 6.2. Let m > a#T/a#. For all i G {0, 1, ■ ■ ■ ,m — 1}, the 

estimate holds 



«#IIK^1ll2,r < C(l|f|li^(0,T;V') + \\9\\hiO,T;Y') + ll^ll 



2,r) 



(37) 



Moreover, if {um}mm is the sequence defined by the step functions Um : / — >■ 
V 



u (t) = I fort = 

then there exists u such that 

Um ^ u in L^(0, T; V). 



Proof. Testing in (36) witli v = and using (11), we get 



a , 



for all 2 G {0, 1, 



m 



a , 



^||Vu^+l||2 



2,Q 



+ 



r 2h 



1}. Hence, applying ^ and Lemma 6.1 it follows 

P^'\\v' + Cy\\g''-'\\Y'y + 



2a# 



^iir„.i+iiii2 
' 2 



with Cy standing for the continuity constant of iJp^(i7i) ^ Y. Summing on 
k = 0, ...,i, multiplying by 2h and applying (10), we find 



a#/^J]||u'=+i||^ + a#||[«^+^]||^,r<-/^E(l|f'llv' + C^yll^?'lly') + 

Cji: 

k=0 ^ 



k=l 



-a#hY,\\[u'^']\\lr + c^*\\S\\lr- 



k=0 



Consequently, by the Gronwall Lemma we get (37) and, for z = m — 1, 
<C(||f|| 

L2(0,T;V') + llfi'llL2(0,T;y') + ll'S'L.r)' (38) 

Thus we can extract a subsequence, still denoted by u^, weakly convergent 
to u e L2(0,T; V). □ 
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Proposition 6.3. Let m > a^T/a^ and Um '■ [0,T[— )■ L'^(T) be defined by 



Umit) = < 



h 



h 



fort = 



on r. 



If the assumptions ^-(11), (I4) and (20)- (22) are fulfilled, then the estimate 
holds 

+ WUmWi^ < C(||f ||i.(o,T;V') + \\9\\Uo,T-x') + ll^llv)- (39) 
Hence, we can extract a subsequence, still denoted by Um-, weakly convergent 

toU e 

Proof. For a fixed t, there exists i G {0, ■ ■ ■, m — 1} such that t G]ti,m; ^i+i,m]- 
Choosing v = u* as a test function in (36), we have 



n Jt Jt 

< {g^+\ u\ - u\^\ + J^ji[u'])ds + (r+\ u^+i - u%. 

Summing on k = 0, i and remarking that 



fc=0 



/n 



In 



In 



dx 



then we find 

fc=0 



0" 



# 



-(f\ u°)o - ^(f^'+i - f^ u'^)n + (r+\ u^+i)^,. (40) 



fc=i 
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Using (|20|)-(|2T}, it follows 

i i 

5^(f'^+l-f^u'=)^<rfl/^5^||u'=||v; 



k=l 



k=l 



J2{9''''-9\n1)r<d2hCYj2\\^'\W- 



k=l 



k=l 



Therefore, inserting the above inequalities in (40), applying (38) and gath- 
ering (37), it results (39). □ 

We again have to relate the weak limits u and U. 
Proposition 6.4. Let u and U he the weak limits obtained in Propositions 



6.S\ and \67!^ respectively. Then 

dt[u] = U in L'^{T), for almost all t E I. 

Proof. For a fixed t, there exists i G {0, ■ ■ ■, m — 1} such that t Gjtj,™.; ti+i,rn\- 
By construction 



Um{r)dT = ^ 



i-l .ik+l)h rk+ll 



k=0 



kh 



h 



-dr + 



ih 



h 



-dr on r. 



Setting the Rothe sequence {u^jmeN defined by 

Um{x,t) = u\x) + {t- ti^rn) : — in L 



h 



for alH G {0, 1, ■ ■ ■ , m — 1} (compare to Definition 3.1 ) under m > a^T/a#, 
it results ^ 

/ Um{T)dT = [Urn]{t) - S OU T. 

Jo 

From the Riesz theorem we get 

{[um]{t) - S,v) = [ {Um{r),v)dT, E L^T) . 
Jo 

Indeed, the right hand side of the above equation is a bounded linear func- 
tional in L'^iT), representable thus (uniquely) by the element [um]it) — S from 
L^(T). Also there exists w G C([0,T]; ^^(r)) such that 

{w{t),v) = [ {U{T),v)dT, yveL^{T). 
Jo 
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lim {[um]it) - S - w{t),v) 



Then we have 

hm [ {UUr)-U{r),v)dT = 0. 

Let us prove that the norms of the functions [um] are uniformly bounded 
with respect to t G / and m. From the estimates (37) independent on i and 
m, and considering 



IM(t)||2,r=||M 1 



t-ti 



2,r 



then, we get 



Ur. 



li°°{0,T;L2(r)) ^ C'(||f ||^2(o,T;V') + II fi' II L2{0,T;y') + ll'S'||2,r) 



Hence, the Lebesgue Dominated Convergence Theorem yields 



lim 

m—^+oo 



i[um]{t) - S -w{t),v)dt = 0, \fv e L'^iT). 



Proceeding as in the proof of Proposition 3.4, we end up with 



[vm-s 



U{T)dT. 



□ 



6.3 Passage to the limit on m — > +oo 

Denoting fm(t) = f'"^^ and g„,{t) = 51*+^ for t e]ti^m,ti+i,m] and z G {0, ■ 
■,m — 1}, we have 



aVum ■ Vvdxdt - 
> / aW\irr,\'^dxdt 



{gm,vi - Umi)rdt + / aUm[v]dsdt+ / j{[v])dsdt> 

aUm[Um]dsdt+ / j{[Um])dsdt+ / {fm,V -Um)Qdt. 

Jt, Jo 



From Propositions |6.2| and |6.3| to pass to the limit the above inequality 
and recalling the weak lower s.c. property for the first term on the right hand 
side of the above inequality, it remains to prove that 

[um] [u] in 
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Taking — u = — + — u first let us prove that 

[um] - [um] -> in 
Since we have < t — ti^rn < h m. ]ti^rn] ^i+i,m] we obtain 

||[tt^](t) - [Mm](t)||2,r = \\Um\\2,r{h - {t- ti^rn)) < || t^m || 2,r- 



Using (39) we derive 



\[Um] - [Mm]||2,S < 



CT 



m 



|f llL2(o,r;V') + ll5'llL2(0,r;Y') + II II v) """^^ ~^ 0" 



Secondly the Rothe sequence {u^} is bounded in L^(0,T; V), and, from 
Prop. 6.4, the functions dt[ bounded in L'^iT,) then, for a subsequence 

still denoted by [um], the strong convergence holds 

[um] [u] in 

Then it results 







Um[u„i\dsdt — > 



U[u\dsdt 



[dtu\ [u\dsdt. 



Therefore we can pass to the limit to obtain the weak formulation (19). 
From the standard technique to prove uniqueness of solution, the solution u 



to (19) with ([8]) is unique. Then the whole sequence {\im} converges weakly 
to u G L2(0,T; V). 



7 Regularity in time 



Proof of Theorem |2.4[ The proof follows the time discretization argu- 



ment as in Theorem |2.3[ considering the existence of the integral inequality 
(36). Testing in (36) for the solutions u*"*"^ and u* with v = (u*+^ + u*)/2. 



summing the consecutive integral inequalities, and dividing by /;., we deduce 



ha\VZ 



i+l|2 



dx + / a{U'+' - U'p'-^'ds < (f 



+ {9-9 
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with U'+^ = {[u'+^] - [u'])/h on T and Z'+^ = {u'+^ - u')/h G V, and 
taking into account the convexity of j. Applying the relation 2{a — b)a = 
+ (a — 6)^ — 6^ to a = f/*+^ and b = U\ and the assumptions (20)-(21), it 
results 



ha\VZ'+'\'dx+ a{U'+yds< a{Uyds + {di + CYd2)h\\Z'+'\\v. 
n Jt Jr 

Notice that the V-norm can be no equivalent to a seminorm. Thus summing 
on A; = 1, ■ ■ ■ , i (z G {1, ■ ■ ■ , m — 1}) we obtain 



k=l 



(41) 



k=l 



with mh = T. 

Let us determine the estimate for the first term on the right hand side of 



the above inequality. Rewrite the integral identity (36) for i = in the form 



[ (7V(u^ - u°) • V(v - u^)dx + / (tVu° • V(v - u^)cix + 
Jn Jn 

+ 1^ «^^^^(H - [u'])ds + ^{j(H) - Ji[u'])}ds > (f 1 - f (0), V - u^)^ + 
+ (f (0), V - u')n - {g' - 9(0), V, - ul)r - {9(0), vi - ul)r, 



for all V G V, and in particular v = u*^. Thus, we apply the assumption (23) 
with V = and divide by h we deduce 

^ah\VZ'\'dx+l^a (j^-j^^ ds < - f(0)||v' + W - 9(0)^') W^'Wv- 



Then, using ([9])-(10), (20)-(21) and taking the Young inequality into account 
for the right hand side, we get 



tt#m2,r< 



2a# 
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Considering h < a# min{l/cr#, 1} we insert the resulting estimate for 
into (41) concluding 



Hence, applying the Gronwall Lemma Um is uniformly estimated in L°°(0; T; L 
and successively is uniformly estimated in L^(0; T; V). Therefore the ex- 



istence of a solution u G C([0,T]; V) in accordance to Theorem 2.4 can be 
done by similar arguments of passage to the limit (cf. Section |4]). 
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